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I.  INTRODUCTION 


This  report  is  a summary  of  research  carried  out  as  part  of  a 
continuing  study  of  the  interaction  of  underwater  sound  with  the  sea 
floor. 

Under  a wide  variety  of  conditions  of  source-receiver  geometry, 
acoustic  frequency,  sound  speed  profile,  and  water  depth,  sound 
propagation  in  the  ocean  can  be  heavily  influenced  by  the  ocean 
bottom.  Although  sea  floor  structure  and  composition  has  itself  been 
a subject  for  study  by  geophysicists  and  seismologists  for  many  years, 
only  comparatively  recently  has  it  been  appreciated  that  sea  floor 
makeup  can  impact  sound  propagation  over  ranges,  frequencies,  and 
geometries  of  concern  to  Naval  applications. 

The  scope  of  these  applications  which  must  be  concerned  with 
bottom  effects  includes  system  performance  prediction,  system  design, 
interpretation  of  acoustical  data,  the  integration  of  geophysical  data 
into  acoustic  modeling,  and  the  design  of  experiments  intended  to  gather 
acoustic  data  in  the  ocean.  This  broad  range  of  concerns  in  which 
acoustic  bottom  interaction  can  play  a significant  role  requires  various 
levels  of  description  of  bottom  interaction  effects. 

Traditionally,  the  nature  of  bottom  interaction  has  been  characterized 
by  a single  quantity,  bottom  loss.  For  many  applications  knowledge  of  this 
quantity  is  indeed  sufficient,  for  example,  in  ray  trace  estimates  of 
propagation  loss  used  in  some  system  performance  prediction  models. 

However,  more  complex  problems  involving  phase  interference  between  various 
multipaths,  such  as  the  functioning  of  arrays,  the  nature  of  Doppler  line 
broadening,  and  a range  of  problems  concerned  with  multipath,  or  mode 
conversion  due  to  range  changing  bathymetry,  all  require  more  detailed 
characterization  of  the  sea  floor  if  these  phenomena  are  to  be  understood 
quantitatively.  Many  of  these  problems  may  require  fairly  detailed 
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descriptions  of  subbottom  sound  speed  and  attenuation  profiles  as  well 
as  location  of  major  interfaces  (reflectors)  and  possibly  of  shear  wave 
parameters.  Such  problems  go  beyond  a simple  regional  characterization 
of  bottom  loss  estimated  in  1/3  octave  bands,  however  useful  such 
information  may  be  for  certain  problems. 

Applied  Research  Laboratories,  The  University  of  Texas  at  Austin 
(ARL:UT),  has  been  conducting  a study  of  the  influence  of  the  ocean  bottom 
on  sound  propagation  characteristics.  In  view  of  the  various  levels 
of  description  of  the  ocean  bottom  required  by  the  intended  applications, 
this  study  has  encompassed  three  primary  research  areas:  (1)  sensitivity 
of  bottom  loss  to  variations  and  uncertainties  in  subbottom  parameters, 

(2)  the  role  of  the  bottom  in  various  range  changing  environment  problems, 
particularly  slope  coupling  and  bottom  roughness,  and  (3)  bottom  inter- 
action effects  such  as  those  involved  in  array  studies,  cw  line  structure, 
and  the  interpretation  of  experimental  acoustic  data. 

From  the  outset  the  primary  goals  of  this  research  program  have  been 
fourfold:  (1)  to  determine  and  provide  guidance  on  the  level  of  detail 

of  subbottom  parameters  required  for  acoustic  applications  (sensitivity 
studies),  (2)  to  determine  which  aspects  of  mode,  or  multipath  conversion 
caused  by  slope  coupling,  and  roughness,  are  predictable  and  exploitable, 

(3)  to  develop  computational  tools  appropriate  to  the  study  of  a wide 
range  of  complex  bottom  interaction  problems,  and  (4)  to  interact  with 
experimental  measurement  programs  via  exercise  planning  and  data  analysis 
and  interpretation. 

During  FY78,  as  reported  here,  research  is  being  carried  out  on 
problems  associated  with  slope  coupling,  the  effects  of  lateral  changes 
in  subbottom  composition,  bottom  roughness  effects,  and  various 
sensitivity  studies.  Much  of  the  work  during  FY78  has  been  to  develop 
computational  tools  to  deal  with  these  problems,  and  much  of  this  work 
is  expected  to  come  to  fruition  in  FY79.  Examples  of  such  work  would  be 
slope  coupling  problems  approached  via  a mode-mode  coupling  study,  and 
the  effects  of  sediment  shear  waves  on  bottom  reflection  loss. 
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The  remainder  of  this  report  is  divided  into  three  major  sections 
dealing  with  range  changing  environment  problems,  bottom  roughness 
studies,  and  bottom  interaction  sensitivity  studies.  Substantial 
documentation,  in  the  form  of  papers  or  reports,  gives  detailed 
expositions  of  the  bulk  of  this  work.  As  a consequence,  this  report 
attempts  only  to  summarize  the  technical  approach  and  to  give  example 
results.  It  should  be  appreciated  that  every  important  research  area 
discussed  here  is  a subject  of  continuing  work  with  results  to  be 
expected  in  the  directions  indicated  here. 
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II.  PROPAGATION  IN  A RANGE  VARIABLE  ENVIRONMENT 

A.  Motivation 

The  purpose  of  this  year's  work  on  propagation  in  a range  dependent 
environment  has  been  to  lay  the  foundation  for  studying  a number  of 
problems  of  current  interest  involving  sound  propagation  in  a range 
variable  ocean  environment.  These  problems  can  be  divided  into  sloping 
bottom  studies  and  studies  of  the  effects  of  lateral  subbottom  range 
variability. 

Specific  problems  to  be  addressed  include:  (1)  determining  the 
depth  dependence  of  average  propagation  loss  in  the  presence  of  slopes, 
with  applications  to  ambient  noise  and  signal-to-noise  problems, 

(2)  partitioning  of  propagated  energy  between  water  and  bottom  paths, 

(3)  arrival  angle  structure  for  receivers  and  arrays  located  near 
continental  slopes,  and  (4)  the  influence  of  slope  subbottom  structure 
on  slope  enhancement,  bottom  interference  effects,  and  array  performance. 
Theoretical  parameter  studies  and  the  development  of  analytical  and 
numerical  computational  techniques  will  emphasize  use  of  a detailed 
description  of  the  subbottom.  These  studies  will  use  data  from  previous 
and  ongoing  ARL : UT  projects. 

The  work  carried  out  during  FY78,  and  summarized  here,  has 
successfully  produced  the  analytical  and  computational  methods  and 
computer  software  necessary  to  attack  this  broad  range  of  problems.  It 
is  recognized,  and  indeed  essential,  that  at  this  early  stage  certain 
problems  may  be  addressed  with  a level  of  detail  that  is  unnecessary  for 
practical  solutions.  Since  one  of  the  objectives  of  this  work  is  to 
determine  the  level  of  detail  that  will  be  required,  not  only  in  describing 
the  sea  floor  acoustically  but  also  in  modeling  propagation,  the  investiga- 
tion should  proceed  from  a general  point  of  view. 

5 

- - - 


In  the  studies  of  lateral  subbottom  range  variability,  emphasis 
will  be  on  determining  the  sensitivity  of  propagation  characteristics 
to  lateral  variations  in  subbottom  structure.  Such  variations  can 
arise  from  nonhorizontal  layering  beneath  the  water-sediment  interface, 
or  from  gradual  lateral  changes  within  a layer  due  to  differential 
deposition  of  varying  grain  sizes. 

In  order  to  attack  the  problem  of  sound  propagation  in  a region 

of  bottom  range  variability,  one  needs  to  select  the  proper  computational 

approach  to  the  problem.  Since  the  main  concern  is  with  bottom  range 

1 2 

variability,  the  coupled  mode  theory  of  Pierce  and  Milder  was  selected 
since  it  allows  a complete  description  of  the  bottom.  Software  has  been 
developed  for  computation  of  acoustic  fields  in  a range  dependent  environ- 
ment using  coupled  mode  theory,  in  particular,  for  computing  mode  coupling 
coefficients  and  radial  equations  in  the  adiabatic  approximation.  Also, 
work  has  been  done  using  analytical  radial  solutions  and  coupling 
coefficients  in  an  effort  to  understand  the  mode  coupling  process  and  to 
understand  the  physical  phenomena  without  the  massive  computational  effort 
required  in  a general  mode  coupling  solution.  Sensitivity  studies 
dealing  with  the  influence  of  bottom  type  on  the  adiabatic  solution  to  the 
radial  equation  were  also  carried  out  this  year  and  the  results  are 
summarized  in  this  report. 

B.  Theoretical  Background 

The  starting  point  for  this  work  is  the  coupled  mode  theory  of 
1 2 

Pierce  and  Milder.  In  this  theory  one  seeks  to  describe  acoustic 
propagation  in  a range  dependent  medium  within  a formalism  in  which  the 
normal  modes  of  propagation  in  the  waveguide  are  coupled  via  the  range 
dependence  of  the  medium. 

The  development  of  the  coupled  mode  theory  begins  with  the  three- 

— ► 

dimensional  wave  equation  for  a point  source  located  at  xq. 


(1) 


v2«Hx) 
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For  a medium  having  azimuthal  symmetry  and  radial  range  dependence,  the 
solution  to  Eq.  (1)  is  expressed  in  partially  separated  form  as 


Mr.z)  = ^ Rn(r)  ‘t’jjU.O 
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(2) 


where  ip(r,z)  is  the  velocity  potential  and  4>n(z,r)  satisfies  the  following 
differential  equation  at  each  range  point  r: 


' 2 

-^-5-  + k2(r,z)  - k2 (r)  <i>  (z,r)  = 0 

L 3z2  n J n 


(3) 


At  each  range  r,  the  <j>n(z,r)  form  an  orthonormal  set  of  functions 
normalized  such  that 

p<J>  (z,r)  4>  (z,r)dz  = 6 , (4) 

n m n,m 
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where  P is  the  density  of  the  medium. 

If  Eq.  (2)  is  inserted  in  Eq.  (1)  and  the  orthonormality  property 

of  Eq.  (4)  is  exploited,  one  may  obtain  the  following  equation  for  the 

radial  function  R (r). 
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Equation  (5)  is  a set  of  coupled  radial  equations  and  the  A (r)  and 

B (r)  are  referred  to  as  coupling  coefficients.  The  coupling  coefficients 
mn 

are  proportional  to  radial  rates  of  change  of  mode  amplitudes,  and  their 
presence  in  Eq.  (5)  gives  rise  to  a redistribution  of  modal  energies  as 
a function  of  range. 


Since  the  numerical  and  computational  techniques  for  calculating 

the  mode  functions  of  Eq.  (3)  are  well  in  hand,  the  main  effort  involved 

in  the  coupled  mode  approach  is  to  calculate  the  coupling  coefficients 

of  Eqs.  (6)  and  (7)  and  to  solve  the  range  equation  (Eq.  (5)).  The 

coupling  coefficients  A and  B are  perhaps  the  most  important 

mn  mn 

functions  in  the  coupled  mode  theory  since  both  the  numerical  procedures 

for  solving  the  range  equation  and  the  simplifying  approximations  to  the 

range  equation  depend  critically  on  the  range  behavior  of  the  B^  and 

A . Therefore  it  is  necessary  to  be  able  to  compute  these  coefficients 
mn 

as  quickly  and  efficiently  as  possible.  Computer  code  is  now  in  existence 
for  computing  the  coupling  coefficients  of  Eq.  (7)  for  a completely  general 
range  variable  propagation  geometry. 
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If  the  terms  involving  mode  coupling  in  Eq.  (5)  are,  in  some  sense, 
small  compared  to  the  other  terms  of  the  equation,  it  makes  sense  to 
employ  the  adiabatic  approximation  to  Eq.  (5).  Within  this  approximation, 
Eq.  (5)  becomes 


-A?  R + - x-  R + k2(r)  R = 
,2m  rdrm  m m 
dr 


^T1  p(*o>  *m(Vr) 


(8) 


Since  the  radial  equations  in  the  adiabatic  approximation  are  many  times 
easier  to  solve  than  the  coupled  radial  equations,  it  is  extremely  useful 
to  be  able  to  estimate  the  magnitude  of  sediment  range  variations  that 
may  be  tolerated  within  the  adiabatic  approximation.  Computer  code  has 
been  written  to  solve  Eq.  (8)  in  3 general  range  changing  environment. 

An  analytic  solution  to  Eq.  (8)  for  an  isovelocity,  wedge  shaped  ocean 
has  also  been  developed  and  some  results  obtained  through  its  use  will 
be  presented. 

C.  Results  for  Lateral  Variations  in  the  Sea  Floor 


The  first  work  done  using  the  numerical  tools  developed  this 
year  deals  with  the  validity  and  sensitivity  of  the  adiabatic  approxima- 
tion as  a function  of  the  range  dependence  of  sediment  geoacoustic 
parameters  (see  Refs.  3 and  4).  This  work  contains  elements  of  both 
the  sloping  bottom  and  lateral  inhomogeneity  problems. 

There  are  several  criteria  for  the  validity  of  the  adiabatic 
approximation;  the  one  selected  for  use  in  this  work  was  derived  by 
Milder  and  involves  a comparison  of  the  coupling  coefficients  B with 
the  eigenvalue  spacing.  Using  this  criterion,  the  adiabatic  approximation 
is  valid  if 


where  X is  the  mode  cycle  distance  given  by 
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X 
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k - *4.1 
m m+1 


Hence,  if  the  mode  coupling  over  a mode  cycle  is  small,  then  adiabatic 
mode  invariance  is  expected. 

The  first  range  dependent  parameter  investigated  was  the  sediment 
sound  speed.  Figure  1 depicts  the  waveguide  geometry  used  in  the 
numerical  calculations  that  were  performed.  The  first  layer  is  a water 
layer  4500  m deep  having  an  arbitrary  sound  speed  profile  (for  the 
calculations  presented  in  this  paper  a mid-Pacific  type  profile  was 
used).  The  sediment  layer  is  300  m deep  and  is  dependent  on  range  and 
depth  as 


c (r,z)  = c (r)  + g(z-H) 

2 o 

This  investigation  of  the  radial  sediment  sound  speed  gradient  falls 
into  the  lateral  subbottom  inhomogeneity  category. 

Figure  2 depicts  the  relationship  between  mode  number  and  local 
radial  sound  speed  derivative  such  that  Eq.  (9)  is  satisfied  for 
£=0.1.  The  frequency  was  taken  to  be  20  Hz  and  the  sediment  was 
modeled  as  either  clay,  silt,  or  sand.  There  are  about  106  discrete 
modes  for  this  frequency  and  layering  geometry.  Several  interesting 
features  are  evident  in  this  figure.  First,  note  the  large  dip  in  the 
curves  at  about  mode  No.  24.  This  dip  marks  the  dividing  line 
between  modes  with  truning  points  in  the  water  and  those  with  turning 
points  in  the  sediment.  For  modes  with  m<24  one  sees  very  little 
influence  from  the  range  variation  of  the  sediment  sound  speed. 

The  main  findings  of  this  portion  of  this  work  may  be  summarized 


as  follows. 


SEDIMENT 

c 2(r,  i),  p2  c2(r,  z)  = c0(r)  + g(z  - H) 


SUBSTRATE 

c3>  P 3 


FIGURE  1 

WAVEGUIDE  LAYER  AND  SOUND  SPEED  STRUCTURE 
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Based  on  Fig.  2 and  a choice  of  e=0.1  it  can  be  said  that  a radial 

-3  -1 

sediment  sound  speed  gradient  of  about  10  sec  may  be  tolerated 

within  the  adiabatic  approximation  for  all  of  the  sediment  types  and 
for  all  mode  numbers. 

In  the  lower  mode  number  portion  (less  than  m-60)  of  Fig.  2 
the  curves  are  in  the  order  sand,  silt,  and  clay;  in  this  low  mode 
region  sand  can  withstand  the  highest  radial  gradients  followed  by 
silt,  and  then  clay.  This  behavior  is  related  to  the  fact  that  the 
characteristic  acoustic  impedance  of  sand  is  higher  than  that  of  silt, 
which  is  higher  than  that  of  clay.  Hence,  for  the  low  mode  numbers  in 
question,  the  sand  bottom  contains  the  lowest  fraction  of  acoustic  energy 
present  and  is  least  affected  by  the  presence  of  a radial  sediment  sound 
speed  gradient.  Silt  is  affected  to  a greater  degree  than  sand  since 
its  characteristic  impedance  is  less,  followed  by  clay  with  the  smallest 
impedance. 

In  long  range  propagation,  it  is  likely  that  only  the  lower  modes  are 
of  importance.  In  a situation  such  as  this  the  magnitude  of  radial  sound 
speed  gradient  that  may  be  handled  within  the  adiabatic  approximation  can 
be  seen  to  vary  over  an  order  of  magnitude  as  the  bottom  type  varies  from 
clay  to  sand. 

D.  Results  for  Sloping  Bottom  Applications 

The  next  parameter  that  was  examined  was  local  bottom  slope.  The 
waveguide  geometry  used  for  calculation  of  the  normal  modes,  coupling 
coefficients,  and  mode  cycle  distances  is  shown  in  Fig.  3.  There  is 
a water  layer  of  variable  depth  given  by  H(r)  overlying  isovelocity 
sediment  and  substrate  layers.  The  calculations  to  be  presented  were 
performed  at  a range  where  the  water  depth  was  4500  m and  the  sediment 
thickness  was  300  m.  Again,  the  water  sound  speed  profile  is  a typical 
mid-Pacific  type. 
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Figure  4 shows  the  relationship  between  mode  number  and  local  bottom 
slope,  dH/dr,  satisfying  Eq.  (9)  for  e=0.1,  with  the  sediment  being 
treated  as  a clay,  a silt,  and  sand.  For  a given  mode  number  m,  values 
of  dH/dr  lying  beneath  a curve  correspond  to  slopes  that  may  be  treated 
adiabatically , while  those  lying  above  a curve  correspond  to  slopes 
that  may  not  be  treated  adiabatically.  Based  on  the  choice  of  e=0.1. 

Fig.  4 reveals  that  most  of  the  modes  propagate  adiabatically  for  slopes 
up  to  about  0.5°. 

Figure  4 summarizes  our  findings  concerning  the  sensitivity  of 
the  adiabatic  approximation  to  bottom  slope  and  sediment  type.  Note 
that  the  curves  tend  to  lie  in  the  order  clay,  silt,  and  sand  with  clay 
being  able  to  tolerate  more  slope  angle  than  sand.  This  ordering 
follows  the  ordering  of  the  characteristic  acoustic  impedances  of  the 
sediments,  with  clay  having  the  smallest  and  sand  the  largest.  This 
behavior  is  not  inconsistent  with  one's  intuition  that  an  acoustically 
hard,  sloping  bottom  would  give  rise  to  more  mode  conversion  than  an 
acoustically  softer,  more  penetrable  bottom.  The  curves  in  Fig.  4,  for 
the  choice  e=0.1,  indicate  that  mode  conversion  effects  begin  to  arise  in 
sand  for  local  bottom  slopes  between  0.057°  and  0.57°,  whereas  for  clay, 
mode  conversion  effects  begin  at  slope  angles  on  the  order  of  0.5°  to 
1.0°.  Hence,  the  mode  conversion  process  due  to  a sloping  bottom  seems  to 
be  fairly  sensitive  to  bottom  type. 

Another  piece  of  work  completed  this  contract  year  concerns  the 
relation  between  ray  theory  and  mode  theory  concepts  in  a sloping  bottom 
region.  The  basic  question  examined  concerns  the  relationship  of  the 
coupling  coefficients  to  the  physically  meaningful  ray  theory  concept 
of  multipath  conversion.  In  other  words,  is  it  the  presence  of  the 
coupling  coefficients  in  Eq.  (5)  that  gives  rise  to  multipath  phenomena 
or  do  they  play  a more  subtle  role? 

To  address  this  problem,  an  analytic  solution  to  the  adiabatic  range 
equation  for  an  isovelocity  wedge  shaped  ocean  was  derived.  Using  this 
analytic  solution,  it  was  possible  to  extend  the  methods  of  Tindle  and 
Guthrie^  to  perform  a ray  trace  entirely  within  the  context  of  mode 
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theory.  Figure  5 depicts  a ray  trace  for  a source  located  at  the 
surface  at  a point  where  the  water  depth  is  500  m.  The  bottom  was 
treated  as  rigid,  with  a 10°  slope.  The  solid  curve  represents  a 
conventional  ray  theory  ray  trace  and  the  circles  represent  a ray  trace 
performed  using  the  adiabatic  radial  solution.  Note  the  excellent 
agreement  even  for  the  large  (10°)  slope  considered.  This  leads  one  to 
conclude  that  the  adiabatic  solution  describes  the  geometric  properties 
of  the  sound  propagation  while  the  coupling  coefficients  are  responsible 
for  redistributing  modal  amplitudes. 
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III.  BOTTOM  ROUGHNESS  STUDIES 

A.  Review  of  Water-Sediment  Roughness  Studies 

The  problem  of  accounting  for  the  effects  of  bottom  roughness  on 
sound  propagation  characteristics  is  a very  difficult  one  that  has 
received  the  attention  of  many  researchers  over  a period  of  years. 

Since  the  inception  of  this  bottom  interaction  program  a small  effort 
has  been  devoted  to  this  specific  problem.  In  particular,  the  goals 
of  this  work  have  been  to  first  assess  the  usefulness  of  existing 
methods  for  accounting  for  bottom  roughness  and  then  to  carry  out 
sensitivity  studies  using  the  most  promising  of  these  methods.  The  goals 
of  the  sensitivity  studies  are:  (1)  to  determine  the  required  level  of 
detail  of  description  of  subbottom  roughness  (How  much  do  we  have  to  know, 
rms  roughness,  correlation  length,  correlation  function,  etc.?),  and 
(2)  to  assess  the  impact  of  bottom  roughness  effects  on  surveillance 
problems . 

Beginning  in  FY77,  efforts  in  this  area  were  based  entirely  on  a 
smoothed  or  average  boundary  condition  approach;^’ ^ such  an  approach 
was  found  by  previous  study  to  offer  the  besl  method  for  accounting  for 
low  frequency  roughness  effects.  Sensitivity  study  results  concerning 

g 

water-sediment  interface  roughness  were  published  during  FY78.  These 
results  suggested  that  water-sediment  interface  roughness  is  probably 
not  important,  at  low  frequencies,  and  for  high  porosity  sediments 
(clays  and  silts)  such  as  those  found  in  deep  ocean  basins  (abyssal 
plains).  During  FY78  it  was  suggested  that,  for  regions  of  thin  sediment 
cover  such  as  the  Northeast  Pacific,  the  most  important  roughness  effects 
might  occur  at  the  sediment-basement  interface.  Work  was  initiated  during 
FY78  to  extend  existing  theory  to  include  the  sediment-basement  interface. 

i 1 
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Basement  Roughness  Studies 


During  FY78  the  work  of  Kuperman  on  the  randomly  rough  two-fluid 
interface  was  extended  to  the  case  of  a fluid-solid  interface.  The 
aim  was  to  calculate  the  specular  reflection  and  transmission  coefficients 
of  the  coherent  component  of  the  sound  field. 


The  geometry  considered  is  shown  in  Fig.  6.  A two-dimensional 
coordinate  system  is  used  with  the  z axis  increasing  downward.  The  fluid 
(medium  1)  and  the  solid  (medium  2)  are  separated  by  the  interface  given 
by  z=ct(x).  The  origin,  z=0,  is  chosen  so  that  the  average  height  of  the 
surface  is  zero. 


The  exact  boundary  conditions  at  a point  on  the  suiface  (x  ,a (x  ))  are 

o o 

written  in  a coordinate  system  (x',z')  in  Fig.  (6)  whose  z axis  is  normal 

to  the  interface  at  that  point.  The  boundary  condition  is  converted  to 

the  original  coordinate  system  by  means  of  a rotation  by  an  angle 

0 = tan  1 ) and  a translation  along  the  z axis  of  a(x  ).  The 

V dx | xQ / O 

resulting  boundary  conditions  are  expanded  about  z=0  to  second  order  in  a. 

The  next  step  involves  the  separation  of  the  fields  into  mean  and 
stochastic  components.  In  terms  of  scalar  and  vector  potentials  this  is 


*1  = *1  + V1  • 


b — Q + v 
p2  2 2 * 


i|>  = *F  + w , 


where  <p ^ , $ are  scalar  potentials;  is  the  vector  potential  in  the 
solid;  v^,  v^»  and  w are  stochastic  fields  (in  the  variable  x)  due  to 
the  rough  interface;  and  4^,  and  V are  the  mean  fields.  The 
boundary  conditions  can  now  be  ensemble  averaged  to  give 
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FIGURE  6 

GEOMETRY  FOR  BASEMENT  ROUGHNESS  THEORY 
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pl$l  - P 2® 2 + 2p2\ 


ay'  _ a^i 

3x  „ 2 
3x 


) = Fl(w’*';  V1,V2,W) 


a®2  s2 

2 ~ a"x  " r + ^2  = F2(*l**2»,|,«  WW)  * 

*1  - *2  " H = F3(w;  V1,V2,W)  * 


where  the  functions  F^,  F^,  and  F^  are  second  order  in  a and  v^,  v^,  and  w 
are  first  order  in  a.  Subtracting  the  mean  field  boundary  conditions  from 
the  full  boundary  conditions  gives,  to  order  a,  the  boundary  conditions 
satisfied  by  the  stochastic  fields: 


Plvl  ~ P2V2  + 2p2 


ft)  ci(VV’)  • 


S a2, 

2 — r^-  - w"  + - — 

ax 


f = G2  (^1 » ^2  *'F) 


H - g3(vv'')  • 


The  fumtions  G^,  G2 , and  G^  are  first  order  in  a and  independent  of  the 
stochastic  fields. 


At  this  point  Eqs.  (10)  through  (12)  state  that  the  correction  due  to 

2 

roughness  in  the  specular  direction  is  of  order  a because  Eqs.  (10)  through 
(12)  reduce  to  the  plane  interface  boundary  conditions  if  F^=F2=F^=0. 

The  procedure  for  calculating  the  reflection  and  transmission 
coefficients  is  as  follows.  First  set  F1=F2=F3=0  and»  assuming  plane 
waves,  solve  for  the  zero  order  (in  a)  fields  at  the  interface.  These 


22 


I 


f 

can  then  be  used  to  calculate  G^,  G^,  and  G^.  Equations  (13)  through 
(15)  are  then  used  to  obtain  the  stochastic  fields  in  terms  of  the  zero 
order  fields.  The  zero  order  fields  and  the  stochastic  fields  are  then 
used  to  calculate  F , F^ , and  F^.  Equations  (10)  through  (12),  after 
subtracting  the  zero  order  fields,  are  then  used  to  calculate  the  second 
order  correction  to  the  mean  fields.  Assuming  a unity  amplitude 
( incoming  wave  will  yield  the  correction  to  the  reflection  and  transmission 

coefficients. 

The  theory  outlined  above  is  complete  and  in  a form  suitable  for 
< numerical  computation.  Work  is  in  progress  to  develop  computer  software 

for  carrying  out  these  computations.  The  results  of  sensitivity  studies 
to  be  carried  out  using  these  methods  will  be  reported  during  the  next 
fiscal  year. 
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IV.  BOTTOM  INTERACTION  SENSITIVITY  STUDIES 

A continuing  aspect  of  the  bottom  interaction  studies  at  ARL:UT  has 
been  the  use  of  various  "measures"  of  bottom  interaction  as  vehicles  to 
test  the  importance  of  one  or  more  parameter  variations  or  uncertainties. 

In  the  beginning  a commonly  used  measure  was  propagation  loss,  and 
work  was  carried  out  to  determine  the  sensitivity  of  propagation  loss  to 
variations  in  bottom  loss.  More  recent  studies  have  explored  the 
acoustical  importance  of  various  sea  floor  parameters  such  as  density, 
sound  speed,  and  absorption  gradient  using  bottom  reflection  loss  as  a 
measure  of  bottom  interaction. 

Sections  following  deal  with  two  additional  elements  in  this  line  of 
investigation:  one,  the  effect  of  sediment  rigidity,  again  uses  bottom 
reflection  loss  as  the  measure;  the  other  concerns  the  effects  of  bottom 
interaction  on  cw  line  spectra.  A third  section  is  concerned  with  the 
relationship  between  two  commonly  used  measures,  bottom  reflection  loss 
per  bounce  and  mode  attenuation. 

A.  The  Effects  of  Sediment  Rigidity 

1.  Background 

In  FY76  a study  was  begun  of  the  sensitivity  of  bottom  reflection 

loss  to  subbottom  parameter  variations.  Initially  the  work  concentrated 

on  the  properties  of  a fluid  sediment  and  was  carried  out  using  a 

9 

sophisticated  bottom  reflection  loss  model  developed  at  ARL:UT  for  this 
purpose.  This  model  treats  the  sediment  as  a fluid.  The  unique  feature 
of  the  model  is  its  use  of  direct  numerical  integration  to  solve  the  wave 
equation.  This  permits  arbitrary  depth  variations  to  be  treated  accurately 
and  without  resort  to  approximation.  As  a research  tool  the  model  has 
proved  to  be  very  flexible  and  convenient  to  use.  Various  sea  floor 
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parameters  such  as  density,  sound  speed  and  absorption  gradients, 

12  13 

and  substrate  rigidity  ’ have  been  successfully  studied. 


In  FY78  the  direction  of  these  studies  turned  toward  the 
inclusion  of  sediment  rigidity  (shear  waves)  in  the  bottom  reflection 
loss  model.  Guided  by  the  success  of  the  ARL:UT  fluid  model,  it  was 
decided  to  develop  a model  based  on  direct  numerical  integration  for  use 
in  studying  the  acoustic  properties  of  solid  sediments.  The  necessary 
theoretical  work  has  been  completed  and  translated  into  an  operating 
computer  program  which  computes  the  bottom  reflection  loss  for  frequencies 
above  about  10  Hz.  This  model  makes  it  possible,  for  the  first  time,  to 
quantitatively  estimate  the  importance  of  sediment  rigidity  within  the 
context  of  an  inhomogeneous  geoacoustic  model  without  making  the 
approximation  of  many  homogeneous  layers.  A second,  more  complex  (and 
computer  time  consuming)  model  is  being  developed  to  extend  the  frequency 
range  down  to  about  3 Hz  and  verify  the  accuracy  of  the  basic  solid  model 
at  its  low  frequency  end. 


Useful  and  realistic  results  have  been  obtained  by  using 

#ave 
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estimates  of  shear  wave  parameters  from  the  work  of  Hamilton  ’ and 


Shirley  and  Hampton, to  constrain  parameter  ranges  in  sensitivity  studies. 
Initial  studies  of  a hypothetical  turbidite  layer  have  been  completed  and 
suggest  that  sediment  shear  waves  may  be  more  important  than  heretofore 
thought. 

2.  Theoretical  Background 

Theoretical  work  was  necessary  in  two  areas.  The  first  was  the 
determination  of  the  appropriate  form  of  the  depth  separated  wave  equations 
to  be  solved.  Major  considerations  were  the  direct  effect  of  parameter 
gradients  and  the  effect  of  coupling  between  shear  and  compressional  waves 
induced  by  the  gradients.  The  second  area  of  theoretical  work  involved 
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developing  a formalism  for  calculating  the  plane  wave  reflection 
coefficient  which  would  allow  the  use  of  direct  numerical  integration. 


a.  Differential  Equations 

The  model  is  shown  schematically  in  Fig.  7.  A two- 
dimensional,  horizontally  stratified  ocean  bottom  is  considered.  It 
is  composed  of  a solid  sediment  overlying  a solid,  homogeneous, 
semiinfinite  substrate.  Sediment  parameter  values  are  permitted  to  vary 
arbitrarily  with  depth.  Solid  lines  indicate  compressional  waves  and 
dashed  lines  indicate  shear  waves.  The  z axis  increases  downward.  H 
is  the  sediment  thickness  and  0 is  the  grazing  angle.  Subscripts  are 
used  to  identify  quantities  in  different  media:  0 refers  to  the  water 

column,  1 to  the  sediment,  and  2 to  the  substrate. 


At  low  frequencies  the  displacement  u in  the  sediment  can  be 
written  in  a generalized  potential  form  due  to  Richards 


u.  - t 7-~  VP,  (x,z,t)  + ir-TTT 


-1  fx(z)  - 1' 


Fx(8) 


The  scale  functions  f^z)  and  F^z)  satisfy  the  differential  equations 


Fi(z)  V1  /pi  VI 

°1<B)1  + 


fl(z)  p’ 

8i(z)  s 7^0  = ^-Gi(z)  * 


where  is  density,  c^  is  sound  speed,  and  v^  is  shear  speed,  all  of 
which  are  permitted  to  vary  arbitrarily,  but  continuously,  with  depth. 
The  scale  functions  have  the  initial  conditions  F^ (0)=f (0)=1 . A prime 
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FIGURE  7 

SCHEMATIC  OF  THE  MODEL 


is  used  to  denote  a derivative  with  respect  to  z.  Decomposing  the  potentials 


and  (J  , assuming  plane  waves,  gives 


= ^(z)  e 


i (kx-<ut ) 


i(kx-uit) 


5.1  = e 


where  u>  is  the  angular  frequency,  k*kQsin0,  ko=a)/cQ,  and  y is  a unit 

vector  along  the  y axis.  Using  the  equations  of  motion,  the  approach  of 
19 

Richards  can  be  applied  and  extended  to  obtain  the  closed  form  coupled 


equations: 


■I + (Gi  - h - + 37  fi  4 - ko 

/fl\fDl  /D1  \ 

= i Vose^j  ^*1  +\d^'  7*1 


4-1  + (*i  " Gl  “ d^)^i  + pT  (°1  ^2  ‘ kocos2^1 

/Fl\fdl  /d  1 \ 

= -ikocos + !)♦{ 


where  dj=l-(b^/o)2) , D^=l-(B^/<^)  , and 


1 = v?[(pi  8l/pl)  + (2  VI  gl/Vl)  + 8I  ■ 81 
B1  - c2i[(pi  Gi/pi)  + (2  CI  Gi/Ci)  + Gi  ' G1 


i-Gi 
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Equations  (16)  and  (17)  are  equivalent  to  the  equations  of  motion  and 
are  valid  at  all  frequencies.  The  direct  effect  of  gradients  appears  in 
the  and  terms  on  the  left  side  of  the  equations.  The  right  side 
of  these  equations  describes  coupling  between  shear  and  compress ional 
waves . 

Estimates  of  the  magnitude  of  b^  and  for  gradients  typical 
of  marine  sediments  show  that  this  coupling  is  negligible  (d=l,  D=l)  for 
frequencies  above  about  3 Hz.  For  higher  frequencies  we  can  set  d^=D^=l 
and  dj=D^=0  to  obtain  the  uncoupled  equations 


+ K 


ri 


(18) 


where 


♦I  + (8i_Gi)^i  + Bi*i  = 0 • (19) 

= (o/c^2  - k2cos20  » an<* 

- k2cos2G 

Equations  (18)  and  (19)  still  contain  the  terms  due  directly  to  physical 
parameter  gradients. 

Estimates  of  g^  and  show  that  the  second  terms  in  Eqs.  (18) 
and  (19)  can  be  neglected  above  about  10  Hz.  Neglecting  these  terms 
results  in  the  simplified  description  given  by 


+ K - 0 


(20) 


= 0 


(21) 


To  complete  the  description  of  the  acoustic  field,  the 
components  of  the  stress  tensor  which  are  continuous  across  an  interface 
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are  required,  i.e.,  the  normal  stresses  p.  , and  shear  stresses  p. 

jzz  jzx 

These  can  be  written  as 


p . = o . (z)  e 

jzz  jzz 


i(kx-wt) 


, . i(kx-wt) 
p.  = o.  (z)  e 
jzx  jzx 

For  the  water  and  substrate  the  depth  functions  have  the  familiar  form 


o = — p u)  <J>  (z) 
ozz  o o 


o =0 
ozx 


°2zz  = -P2(a)2-2v2k2)  ^2  + 2ip2v2k^2^Z^ 


°2zx  = 2*P2v2k  <t>2^Z^  + P2^u,2-2v2k2^2^Z^ 


The  presence  of  scale  factors  in  results  in  a more  complicated  result 
in  the  sediment  layer.  For  the  high  frequency  uncoupled  limit  of 
Eqs.  (18)  and  (19)  we  have 


j j 2 M 2 

-p1(o>  -2vjk  ) -jr  ~ Pj^GjcJ  -p 


, <1*,  2 

iPjGjcJ  k y-  + 2iPlV;L  k j L- 


2 '•’i  2 T1 

°lzx  ■ -‘"A’l  k 77  + 2iPl  Vlk  77 


. , 2 , 22.  „ 2 1 

+ P1(u  -2v1k  ) — + p181v1  f~ 


These  differ  from  the  usual  results  in  two  ways.  First,  the  potentials 
and  their  derivatives  are  scaled  by  their  respective  scale  factors. 
Second,  new  terms  appear  which  depend  on  the  gradients  of  sediment 
parameters  through  g^  and  G^. 

b.  Mathematical  Formalism 


Most  models  which  include  both  sediment  rigidity  and  the 

depth  dependence  of  sediment  parameters  have  been  based  on  the  Thomson- 
18  19 

Haskell  ’ matrix  method.  In  this  method  gradients  are  approximated 
by  many  thin  homogeneous  layers.  Within  each  layer  the  wave  equations 
are  solved  analytically  in  terms  of  upgoing  and  downgoing  plane  waves. 
The  matrix  approach  allows  the  elimination  of  most  plane  wave  amplitudes 
as  quantities  to  be  calculated,  and  compactly  represents  a structured 
ocean  bottom  in  terms  of  only  four  variables,  one  of  which  is  the 
reflection  coefficient  R. 


Our  approach  is  motivated  by  the  success  of  direct 

numerical  integration  in  modeling  fluid  sediments,  and  by  the  compactness 

of  the  Thomson-Haskell  matrix  representation  of  layered  solids.  Both  of 
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these  can  be  combined  using  the  "propagator"  formalism.  In  this 
formalism  the  differential  equations  describing  shear  and  compressional 
wave  propagation  within  a layer  are  replaced  by  a matrix  relating  boundary 
values  at  the  top  and  bottom  of  the  layer.  The  elements  of  the  propagator 
matrix  can  be  calculated  by  direct  numerical  integration  of  the  depth 
separated  wave  equations. 

The  two  second-order  differential  equations  for  and 
can  be  written  as  four  first-order  differential  equations  for  4> » 4> ^ » 'l'^* 
and  These  four  equations  can  be  written  in  matrix  form  as 


dz 


Ax(z)  *1 


* 
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where  4>^=(4>^,<t>|,<|^,<Jj|)  is  the  four-component  potential  vector  and 
the  matrix  A^(z)  representing  Eqs.  (18)  and  (19)  is 


Since  the  elements  of  A are  continuous  functions  of  depth  the  solution 
20  1 

can  be  written  as 

$1(H)  = E^H.O)  4^(0) 

where  E^(H,0)  is  the  propagator  relating  4^  at  z=H  to  its  value  at  z*0. 
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The  propagator  is  evaluated  by  choosing  a set  of  different 
potential  vectors  at  the  top  of  the  sediment  and  integrating  to  the  bottom. 
A set  of  equations  in  the  unknown  matrix  elements  is  generated  and  solved. 
For  the  differential  equation  sets  of  Eqs.  (18)  and  (19)  and  Eqs.  (20) 
and  (21)  there  are  only  eight  nonzero  matrix  elements  since  the  shear  and 
compressional  fields  are  decoupled.  In  these  cases  only  two  properly 
chosen  Initial  vectors  are  required. 

The  procedure  for  calculating  the  reflection  coefficient 
is  the  same  as  that  of  the  traditional  Thomson-Haskell  matrix  method 
with  the  exception  of  using  the  propagator  to  solve  the  differential 
equations.  Briefly  the  procedure  is  as  follows.  Starting  with  the 
potential  vector  in  the  water  due  to  the  incoming  and  reflected  plane 
waves,  convert  to  the  stress  displacement  vector  and  use  the  continuity 
of  normal  displacement  and  stress  and  shear  stress  to  cross  into  the 
sediment.  Converting  back  to  potentials,  the  propagator  is  used  to  cross 
the  sediment  layer.  Next,  the  sediment-substrate  interface  is  crossed  by 
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converting  to  stress  displacement,  using  the  continuity  conditions  and 
converting  back  to  the  potentials  due  to  the  shear  and  compressional 
plane  waves  in  the  substrate.  This  procedure  results  in  four  equations 
in  four  unknowns,  one  of  which  is  the  reflection  coefficient.  Solving 
the  system  of  equations  leads  to  the  reflection  coefficient.  All  possible 
configurations  of  fluid  or  solid  sediment  and  substrate  have  been  treated 
in  the  matrix  formalism. 


3.  Computation  of  Bottom  Loss 

A model  based  on  the  inhomogeneous  wave  equations  (Eqs.  (20) 

and  (21))  has  been  developed  into  a functioning  computer  program.  This 

model  is  the  most  efficient  in  terms  of  computation  time,  and  provides 

accurate  results  for  frequencies  above  10  Hz.  The  numerical  integration 

9 

technique  has  been  discussed  in  detail  elsewhere. 

Efforts  are  currently  under  way  to  implement  a model  based  on 
Eqs.  (18)  and  (19).  This  will  be  used  to  check  the  accuracy  of  the  basic 
model  at  low  frequencies,  and  also  to  extend  the  model  capability  to 
yet  lower  frequencies. 


In  the  basic  model  at  a given  frequency  and  grazing  angle  the 
reflection  coefficient  is  calculated  for  the  appropriate  case  of  fluid 
or  solid  sediment  or  substrate.  The  quantity  delivered  as  output  is  not 
R itself  but  rather  the  reflection  loss  (bottom  loss). 


RL  = -20  log10(|R|) 


and  the  phase  angle  of  R,  6=tan  ^ [ Im(R) /Re (R) ] . These  are  available 
both  as  hard  copy  printout  and  as  a disk  file  prepared  for  plotting. 

Two  types  of  output  are  possible:  (1)  RL  as  a function  of  grazing  angle 
at  a given  frequency,  or  (2)  RL  as  a function  of  frequency  at  a given 

grazing  angle.  In  this  latter  case  the  attenuation  is  assumed  to  be 

21,22 

* ' 


proportional  to  frequency 


and  is  input  normalized  to  the  value  at 


34 


1 kHz.  Because  of  a local  error  tolerance  feature,  execution  times 
increase  with  wave  frequency  as  wavelengths  become  small. 

Typical  input  parameters  are  shown  in  Table  I for  the  sediment 

23 

structure  of  Fig.  8 taken  from  Fryer.  The  depth  grid  on  which  the 
parameters  are  specified  is  arbitrary.  Linear  gradients  are  assumed 
between  grid  points.  No  assumptions  are  made  about  the  relationship 
between  parameters. 


Figure  9 shows  typical  results  obtained  for  the  profile  of 

Fig.  8.  Three  cases  at  a frequency  of  20  Hz  are  shown.  For  a fluid 

substrate  no  difference  was  seen  between  a solid  (dashed  curve)  and  fluid 

(dotted  curve)  sediment.  For  a fluid  sediment  and  solid  substrate  the 

expected  increase  in  RL  is  seen  between  0 =50°  and  0 =70°,  where  0 and  0 

s p ’ s 

are  the  critical  angles  for  shear  waves  and  compressional  waves  in  the 
substrate.  For  a solid  sediment  and  a solid  substrate  (solid  curve)  the 
resonance  structure  between  0 and  0 is  modified.  No  other  major  effect 


P 


The  correct  operation  of  the  program  was  verified  as  follows. 

For  a fluid  sediment  the  results  for  solid  and  fluid  substrate  conditions 

9 

were  compared  with  those  obtained  from  an  existing  model,  which  treats 
fluid  sediments  only.  For  a solid  sediment  the  results  for  fluid  and 
solid  substrates  were  compared  with  those  obtained  for  particular  choices 
of  depth  dependence  for  which  analytic  solutions  of  the  wave  equations 
are  possible.  The  reflection  coefficient  was  calculated  by  solving  the 
set  of  six  (fluid  substrate)  or  seven  (solid  substrate)  simultaneous 
equations  generated  by  the  continuity  conditions  for  a homogeneous  sediment. 
The  comparisons  showed  excellent  agreement  and  thus  verified  the  logic  flow 
of  the  program  and  the  accuracy  of  the  numerical  integrator.  In  the  course 
of  operation  one  area  of  inaccuracy  has  been  noted.  For  small  grazing 
angles  and  deep  sediments,  a loss  of  precision  is  possible  because  of  the 
exponential  nature  of  the  wave  functions.  However,  this  is  not  a serious 
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TABLE  I 

PHYSICAL  PARAMETERS  OF  SEDIMENT  STUDIED.  THE  ATTENUATIONS 
ARE  GIVEN  FOR  A FREQUENCY  OF  20  Hz. 


Depth 

m 

c 

o 

m/  sec 

V 

o 

m/sec 

P 

g/ cm^ 

a 

P 

dB/m 

J 

a 

s 

dB/m 

Water 

1530 

i 

i 

i 

1.03 

1 

1 

1 

i 

i 

i 

0 

1510 

116 

1.53 

0.0013 

0.169 

36 

1582 

283 

1.579 

0.0020 

0.112 

120 

1674 

391 

1.689 

0.0040 

0.172 

518 

1992 

621 

2.010 

0.0027 

0.087 

Substrate 

4460 

2400 

2.46 

0.00016 

0.00079 

COMPRESSIONAL  WAVE  SPEED  - m/sec 
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problem,  since  the  wave  does  not  penetrate  deeply  into  the  sediment  and 
the  thickness  of  the  layer  can  be  reduced  until  accurate  results  are 
obtained. 

4.  Results 

24 

This  model  has  been  used  by  Vidmar  and  Foreman  to  study  the 

effect  of  sediment  rigidity  on  R for  the  hypothetical  turbidite  layer 

23 

of  Fig.  8 which  was  obtained  from  Fryer.  These  parameters  represent 
realistic  depth  variations  as  deduced  from  geophysical  data.  There  are 
also  depth  dependent  attenuations  not  shown. 

Figure  9 shows  a comparison  of  the  calculated  reflection 
coefficients  of  this  sediment  for  three  different  bottom  configurations 
at  a frequency  of  20  Hz.  The  dotted  curve  is  obtained  by  setting  the 
shear  parameters  to  zero  and  corresponds  to  the  case  of  a fluid  sediment 
and  a fluid  substrate.  The  oscillations  are  due  to  the  interference 
between  the  wave  reflected  directly  from  the  water-sediment  interface 
and  the  wave  penetrating  the  sediment  and  returning  to  the  water.  To 
isolate  the  effect  of  substrate  rigidity,  the  dashed  curve  was  generated 
for  the  case  of  a fluid  sediment  but  a solid  substrate.  Comparison  with 
the  dotted  curve  shows  that  substrate  rigidity  is  responsible  for  the 
increased  loss  and  the  structure  between  50°  and  70°.  The  critical 
angle  for  shear  wave  propagation  in  the  substrate  is  50°,  and  70°  is 
the  critical  angle  for  compressional  waves.  The  increased  reflection 
loss  between  the  two  critical  angles  is  due  to  energy  carried  away  by 
the  substrate  shear  waves.  The  effect  of  sediment  rigidity  can  now  be 
seen  by  comparing  the  solid  and  dashed  curves.  The  solid  curve  was 
generated  for  a solid  sediment  and  substrate.  The  difference  between 
these  two  curves  is  entirely  due  to  sediment  rigidity.  There  is  a 
negligible  effect  except  between  the  two  critical  angles;  here  sediment 
rigidity  produces  a shift  in  the  location  of  the  peaks. 
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| Figure  10  shows  the  effect  of  sediment  rigidity  when  the 

sediment  thickness  is  reduced  to  252  m,  still  at  20  Hz.  Both  curves 
were  obtained  for  a solid  substrate.  As  in  the  previous  and  following 
figures,  the  solid  curve  is  for  solid  sediment  and  the  dashed  curve  for 
fluid  sediment.  In  addition  to  the  expected  peak  shifts  between  the 

i 

critical  angles,  sediment  rigidity  has  the  additional  effect  of 
1 increasing  the  reflection  loss  above  24°. 

I 

As  the  sediment  thickness  is  further  reduced  a pattern 

develops.  There  is  a minimum  angle,  0q,  above  which  sediment  rigidity 

increases  the  reflection  loss.  The  angle  0q  steadily  decreases  as  the 

thickness  decreases.  The  situation  at  36  m thickness  is  shown  in  Fig.  11, 

still  at  20  Hz.  Here  0 has  decreased  to  zero.  Sediment  rigidity  clearly 

o 

dominates  the  reflection  loss  at  low  angles  with  a dramatic  increase  of 
20  dB  near  13°. 

The  effect  of  sediment  rigidity  is  expected  to  decrease  at 
higher  frequencies.  Figure  12  shows  the  reflection  loss  at  200  Hz  for 
the  same  36  m layer.  The  increased  reflection  loss  is  indeed  smaller, 
but  is  not  yet  negligible.  Even  at  this  relatively  high  frequency  there 
is  an  increase  of  about  4 dB  between  15°  and  30°. 

Several  observations  can  be  made  based  on  these  results. 

First,  sediment  rigidity  does  matter.  Second,  the  reflection  loss  depends 
on  layer  thickness  in  a regular  fashion.  The  angle  above  which  sediment 
rigidity  is  important  increases  as  the  thickness  increases.  Thirdly, 
excitation  of  shear  waves  at  the  water-sediment  interface  is  not  important, 
whereas  conversion  of  compressional  waves  at  the  substrate  interface  is 
dominant . 


This  last  point  requires  some  elaboration.  Significant  changes 
in  the  reflection  loss  occurred  when  only  the  sediment  thickness  was 
changed.  Since  the  upper  interface  parameters  were  constant,  physical 
processes  occurring  there  are  not  likely  to  be  the  major  cause  of  the 


40 


J 


I 

I 


i 

i 

i 


gp  - sson  Nouoauaa 


USE 

=>  O P 

“*  g* 

• y 


UJ  H 
0£ 


SOLID  SEDIMENT 
FLUID  SEDIMENT 


9P 


ARL:UT 

AS-78-1739-S 

PJV-GA 

11-27-78 


! 

{ 

I 

I 


; 


increased  reflection  loss.  If  we  consider  the  interaction  of  the 
compressional  wave  with  the  substrate,  we  notice  that  there  is  an  angle 
above  which  the  interaction  is  strong.  For  thick  sediments  the  angle  is 
large.  As  thickness  decreases,  the  angle  decreases.  This  is  qualita- 
tively the  same  dependence  found  for  0q.  This  similarity  suggests  the 
importance  of  compressional  wave  conversion  at  the  lower  interface. 
Preliminary  work  on  a small  parameter  expansion  of  the  reflection  and 
transmission  coefficients  at  the  interfaces  also  verifies  this  conclusion. 

An  additional  application  of  this  work  during  FY78  was  the 
estimation  of  the  importance  of  shear  wave  effects  to  the  development 
of  a geoacoustic  model  for  the  BEARING  STAKE  area.  This  work  was 
carried  out  in  support  of  the  BEARING  STAKE  exercise  data  analysis. 

B.  Bottom  Interacting  Multipath  Effects  on  Line  Structure 

The  presence  of  a significant  spread  in  multipaths,  and  hence 
group  velocities,  will  cause  a moving  cw  source  to  produce  Doppler  effects 
which  can  be  significantly  different  for  the  various  multipaths.  As  a 
result,  the  spectra  associated  with  such  acoustic  fields  will  display 
Doppler  broadening  effects.  Although  such  an  effect  will  occur  whenever 
more  than  one  multipath  is  significant,  an  important  spread  in  multipath 
group  velocities  will  usually  occur  only  when  a number  of  bottom  inter- 
acting paths  are  present.  Thus,  the  structure  of  cw  line  spectra  provides 
an  additional  measure  of  the  importance  of  bottom  interaction,  and  one 
which  is  of  considerable  practical  importance. 

25 

A general  purpose  method  was  recently  developed  at  ARL:UT 
to  deal  with  the  problem  of  Doppler  broadening  and  the  interplay  of 
propagation  effects  and  signal  processing  parameters.  This  method  relies 
on  a normal  mode  theory  of  the  field  generated  by  a moving  source  and 
includes  a complete  structured  bottom  as  well  as  arbitrary  source-receiver 
geometry.  The  remainder  of  this  section  will  be  devoted  to  exhibiting  a 
number  of  results  obtained  using  this  model  in  an  effort  to  show  that  the 
anticipated  sensor  geometry,  source  speed,  and  integration  time  dependencies 
are  indeed  obtained  quantitatively. 
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1.  Receiver  and  Source  Depth  Effects 


Consider  now,  as  an  example,  a source  moving  on  a constant 
course  at  constant  speed  v,  passing,  at  CPA,  within  30  km  of  a receiver. 

The  geometry  is  shown  in  Fig.  13  where  the  source-receiver  range  R(t)  will 
be  given  by  (30  km)  x cos9(t).  Figures  14  and  15  show  the  results  of 
predicted  25  Hz  line  spectra  for  a 5 m source  depth  and  1000  m and  4500  m 
receiver  depths,  respectively.  The  sound  speed  profile  used  in  these  calcu- 
lations is  taken  from  the  CHURCH  OPAL  exercise  area  with  a bottom  descrip- 
tion appropriate  to  this  area  of  the  North  Pacific.  The  spectra  shown 
here  have  1 mHz  resolution  (1000  sec  integration  time).  The  multipath 
broadening  effects  are  quite  evident,  as  is  the  basic  Doppler  shift 
(upward  from  25  Hz).  Because  the  various  multipaths  interfere  in  different 
ways  at  the  two  receiver  depths  and  because  at  the  deep  receiver  (380  m 
above  the  bottom)  there  will  be  a higher  proportion  of  bottom  interacting 
energy,  the  details  of  line  structure  are  different  at  these  locations. 


Figures  16  and  17  show  the  same  configuration  except  that  the 
source  depth  is  100  m rather  than  5 m.  Again  it  is  observed  that  receiver 
depth  affects  line  structure.  Moreover,  comparing  Figs.  11  and  12,  we  see  that 
the  low  Doppler  bottom  interacting  modes  are  more  important  for  the  shallow 
source  than  for  the  100  m source. 

2 . Range  Effects 

Figures  18  through  21  show  a sequence  of  spectra  for  a 100  m source 
for  ranges  of  200,  100,  50,  and  10  km  from  CPA  (again,  30  km).  It  will  be 
observed  that  as  the  source  moves  closer  to  the  receiver  the  bottom 
interacting  modes  become  mere  important,  leading  to  more  pronounced 
Doppler  shifts  for  these  paths.  The  receiver  depth  in  these  four  cases  is 
4833  m (50  m above  the  bottom).  It  is  clear  from  these  calculations  that, 
especially  for  cases  of  receiver  near  a boundary,  line  shape  will  also 
depend  on  bottom  loss 
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FIGURE  13 

SCHEMATIC  ILLUSTRATION  OF  SOURCE  TRACK 
AND  RECEIVER  GEOMETRY 
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FIGURE  14 

Normalized  Spectrum  at  a 1000  m Receiver  for  a 25  Hz,  5 m Deep 
Source  Located  20  km  Downrange  from  CPA 


Receiver  Depth:  4500  m 
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FIGURE  15 

Normalized  Spectrum  at  a 4500  m Receiver  for  a 25  Hz,  5 m Deep 
Source  Located  20  km  Downrange  from  CPA 


FIGURE  16 

Normalized  Spectrum  at  a 1000  m Receiver  for  a 25  I 
Source  Located  20  km  Downrange  from  CPA 


Receiver  Depth:  4500  m 
Source  Depth:  100  m 


Receiver  Depth:  4833  m 
Source  Depth:  100  m 
Bottom  Depth:  4883  m 


FIGURE  19 

Spectrum  at  a Receiver  50  m Above  the  Bottom 
for  a 25  Hz  Source,  100  km  Downrange  from  CPA 


Receiver  Depth:  4833  m 
Source  Depth:  100  m 


FIGURE  20 

Spectrum  at  a Receiver  50  m Above  the  Bottom 
for  a 25  Hz  Source,  50  km  Downrange  from  CPA 


Receiver  Depth:  4833  m 
Source  Depth:  100  m 


FIGURE  21 

Spectrum  at  a Receiver  50  m Above  the  Bottom 
for  a 25  Hz  Source,  10  km  Downrange  from  CPA 


3.  Source  Speed  Effects 


Figure  22  shows  the  same  situation  as  in  Fig.  15  except  that 
the  source  speed  has  been  changed  to  15  kt.  It  will  be  observed  that  not 
only  is  the  overall  Doppler  shift  increased  (the  rms  frequency  is  changed 
from  25.020  to  25.053  Hz)  as  would  be  anticipated  on  the  basis  of  elementary 
Doppler  shift  theory,  but  the  line  structure  itself  is  altered.  However, 
both  of  the  rms  Doppler  shifts,  20  mHz  and  53  mHz,  and  the  rms  line  widths, 
3.2  mHz  and  9.6  mHz,  are  increased  by  approximately  a factor  of  3 when  the 
source  speed  is  increased  from  5 to  15  kt. 

4.  Frequency  Effects 

It  can  be  expected  that  there  will  be  some  changes  in  line 
structure  between  lines  of  appreciably  different  frequencies,  particularly 
when  bottom  interacting  multipaths  (and  therefore  frequency  dependent  bottom 
loss)  are  important.  Figure  23  shows  a situation  identical  to  that  of 
Fig.  16  except  that  the  line  frequency  has  been  changed  from  25  to  100  Hz. 

It  will  be  observed  that  the  secondary  lines  seen  in  Fig.  16  are  absent 
in  Fig.  23.  Additional  possibilities  for  frequency  effects  lie  in  the 
areas  of  amplitude  and  frequency  fluctuations  of  the  spectrum. 

5.  Integration  Time  Effects 

It  would  of  course  be  expected  that  decreased  integration 
time  (decreased  resolution)  would  act  to  broaden  lines  and  "smear  out" 
secondary  lines.  Direct  calculations,  of  the  type  shown  in  this  section, 
have  been  carried  out  and  verify  this  expectation.  However,  there  is 
another  integration  time  related  effect  which  is  also  contained  in  the 
model  predictions,  and  which  acts  in  the  opposite  direction.  When  the 
rate  of  change  of  range  is  large,  the  rms  Doppler  shifted  frequency  also 
changes  rapidly  with  time.  This  results  in  an  FM-slide  type  of  effect, 
whatever  the  multipath  structure.  For  sources  moving  at  fairly  high 
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Receiver  Depth:  1000  m 
Source  Depth:  100  m 

Bottom  Depth:  4883  m 


FIGURE  22 

Normalized  Spectrum  for  a 15  kt  Source 
20  KM  FROM  CPA,  AT  A 1000  M RECEIVER 


Receiver  Depth:  4500  m 


AS-79-372 


FIGURE  23 

Normalized  Spectrum  for  a 100  Hz  Source 
20  KM  FROM  CPA,  AT  A 4500  M RECEIVER 


speeds,  this  effect  can  result  in  very  broad  lines  with  a shape 
completely  characteristic  of  an  FM  slide. 

Under  these  conditions,  a decrease  in  integration  time  will 

result  in  a narrowing  of  the  line,  not  a broadening.  Calculations 

have  been  carried  out  confirming  the  existence  of  this  effect  for  cases 

of  practical  importance.  Thus,  line  width  is  dependent  on  (1)  integration 

time,  (2)  multipath  structure  of  the  field,  and  (3)  the  range-rate  effect 

near  CPA.  Since  the  combined  effects  of  (1)  and  (3)  can  act  in  a way 

that  is  opposite  to  the  combined  effects  of  (1)  and  (2),  the  factors 

determining  line  widths  in  a particular  case  can  be  quite  complex.  These 
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issues  have  been  discussed  in  a recent  paper  by  Hawker  and  Shooter. 

C.  The  Relationship  Between  Mode  Attenuation  and  Bottom  Loss/Bounce 

In  the  context  of  propagation  problems  the  loss  effects  of  the 
ocean  bottom  are  taken  into  account  in  two  principal  ways:  (1)  mode 
attenuation  in  normal  mode  models,  and  (2)  bottom  loss  used  in  ray 
models.  Of  course,  bottom  loss  could  be  computed  from  geoacoustic 
description  using  wave  theory  methods,  as  is  done  in  reflection  loss 
models,  and  the  result  used  in  a ray  theory  propagation  model. 

Nevertheless,  the  disjoint  approach  to  loss  effects  brought  on  by 
a reliance  on  either  bottom  loss  per  bounce  or  mode  attenuation  is  not 
only  conceptually  troublesome,  but  it  has  inhibited  progress  in  bottom 
loss  measurement  and  in  parameter  studies  of  bottom  interaction  effects. 

During  the  present  year  an  investigation  of  this  subject  was 

carried  out  with  the  purpose  of  obtaining  relationships  between  mode 

attenuation  (and  normal  mode  quantities  such  as  mode  interference 

distance)  and  bottom  loss  per  bounce.  The  principal  results  of  this 

27 

study  are  given  by  Tindle  with  additional  theoretical  background  by 
28 

Weston  and  Tindle. 
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The  primary  results  of  this  study  are  twofold: 

. The  mode  attenuation  per  mode  interference  distance 

D =2 u (dn/dk  ) given  by  6 / D is  essentially  identical  with 
n n n n 

reflection  loss.  This  equivalence  has  been  shown  within  the 
context  of  WKB  theory  as  well  as  generally,  using  numerical 
computations. 

. The  mode  interference  distance  was  also  found  to  be  given  by 

4ylkn 

D°"(n  vH>)2+k<H>)2 


where  U and  U'  are  the  values  of  the  normal  mode  evaluated 

n n 2 2 2 
at  the  bottom  of  the  water  column,  and  y^=k  -kn« 

These  results  are  important  for  several  reasons  going  beyond 

the  conceptual  advances  they  produce.  The  combination  of  these  two 

results,  together  with  the  identification  of  U'(H)/U  (H)  with  a particular 

n n 

function  of  the  reflection  coefficient,  opens  the  possibility  of  normal 
mode  computation  based  on  bottom  loss  rather  than  a geoacoustic  subbottom 
description.  Although  a major  thrust  of  this  overall  bottom  interaction 
research  program  has  been  to  emphasize  the  utility  of  a geoacoustic 
description,  the  existence  of  a large  body  of  measured  bottom  loss  data 
requires  that  bottom  loss  be  a possible  input  in  propagation  models. 
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V.  SUPPORT  ACTIVITIES 

With  the  analytic  intractability  of  many  of  the  acoustic  equations, 
computer  assisted  numerical  computation  serves  a vital  function  in  the 
bottom  interaction  program.  The  acquisition  of  the  CYBER  171  by  ARL:UT 
in  January  1978  has  greatly  improved  the  computing  capabilities  of  this 
laboratory.  Tasks  formerly  regarded  as  infeasible  due  to  limited  core 
and  disk  space  and  excessive  computation  time  can  now  be  performed. 

During  1978,  the  major  software  tools  previously  developed  by  the 
bottom  interaction  program,  including  the  ray  and  mode  models  and  the 
bottom  loss  reflection  model,  were  converted  to  CYBER  compatible  FORTRAN 
Some  of  these  models  were  restructured  and/or  expanded  to  enable  more 
efficient  use  of  the  computer. 

The  major  software  development  during  this  year  was  an  extension 
of  the  existing  BOTLOSS  (bottom  loss  reflection)  model.  Previously 
dealing  only  with  fluid  sediments,  this  model  will  now  compute  effects 
due  to  shear  waves  in  the  sediments.  The  program  has  been  extensively 
tested  and  is  being  used  in  production  runs.  The  software  package 
centered  around  the  normal  mode  model  was  enlarged  by  the  design  and 
implementation  of  a model  calculating  the  time  series  and  line  structure 
generated  by  a moving  source. 

A laboratory-wide  computer  terminal  network  was  set  up  in  April, 
and  with  this  means  of  easy  access  to  the  CYBER,  usage  of  all  software 
was  greatly  increased.  It  became  apparent  that  most  users  were  unable 
to  take  advantage  of  the  sophisticated  models  due  to  their  incomplete 
knowledge  of  programming,  job  operations,  and/or  specific  model  require- 
ments. To  remedy  this  and  to  provide  easy  and  effective  use  of  all 
production  models,  a model  operating  system  was  initiated.  This  system 
facilitates  the  running  of  the  models,  reduces  the  possibility  of  errors, 
and  requires  only  a very  limited  knowledge  of  CYBER  job  control.  All 
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bottom  interaction  software  has  been  incorporated  into  this  system,  and 
consequently  has  been  used  extensively.  Computer  accessible  documentation 
directing  the  use  of  this  system  has  been  installed.  A library  of  sound 
speed  profiles  is  also  available  via  terminal  interaction. 
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